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Abstract 


A new, geometrical, first order, nonresonant, frozen orbit theory has 
been developed based on Orlov's uniformly rotating plane of constant 
inclination. Perturbation spectra generated from a 90th order subset of 
OSU86F are shown for the ill-fated 1984 JHU/APL SAGE proposal for a pair of 
TRANSIT satellites at 400km altitude with a 93!5 inclination. 


Introduction 


To perform the integration for a geodetic satellite a new form has 
been developed for the geopotential on an orbit as 
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where the symbols are defined in [Melvin, 1988a] and an algorithm is 
derived for the generation of the complex coefficients from a spherical 
harmonic geopotential model. 

The geopotential variation 
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is the distance of an equipotential surface from sphere of radius,^ 
spectrum in Fig. 1 is a log- log bar graph of the amplitudes 2a |R 
the frequency is |pn + q$| . The along- track deflection evaluated 
nominal circular orbit is 
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and its spectrum is plotted in Fig. 2. 
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Position Perturbations 

The real forms of the series are used for the integration in [Melvin 
1988a] but it is apparent from [Melvin, 1988b] that the formulas are much 
simpler with complex coefficients. In fact the orbit perturbations are of 
the same form as the disturbing potential 
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where the coefficients are computed from the algorithm 
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After multiplication by 2a, the amplitudes of the coefficients of (3), (5) 
and (7) are plotted in Figs. 3, 5, and 7 from which it is seen that 
position perturbations accentuate the low and near orbital frequencies and 
attenuate the high frequencies. Expressions for the cross-track (3), 
radial (5), and along-track (7) position perturbations from the Kaula-Allan 
theory are found in [Rosborough and Tapley, 1987]. 


Velocity Perturbations 

By use of coordinates at the nominal satellite position, the cross- 
track, radial, and along-track velocity components are 

v = a(<p - ftsintcosu) - a[ £ X i(prt+q&)<p e*^ U “ ftsimcosu] , (4) 
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The amplitudes of the coefficients of (4), (6) and (8) are plotted in Figs. 
4, 6, and 8. By comparison of the high frequency fall-off in the position 
and velocity spectra, it is clear why more geopotential information is 
obtained from Doppler beacon satellites than from skin tracking even if it 
is by laser reflectors . 


Intersatellite Measurements 

For a close satellite pair in same orbit, the intersatellite range is 
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A time derivative yields the intersatellite range rate as 
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For a nominal separation of a<5u=100km the spectra of (9) and (10) are 
plotted in Figs. 9 and 10. A comparison of the high frequency portions of 
Figs. 2 and 10 and the foregoing formulas show 
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from which the result of [Comfort, 1973] is modified to state that at high 
frequencies for a pair of satellites flying in formation the intersatellite 
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range rate mimics the along track deflection. The paucity 
lines near twice per orbit in Fig. 10 indicates that complete 
recovery is not possible with only intersatellite range rate 


of spectral 
geopotential 


Gravity Gradient 

. E ^ ua Py simple formulas for the gravity gradient tensor along a 
nominal circular orbit in Orlov's plane are found in [Melvin, 1988b] as 
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Geopotential Recovery 

^ iS . sub j ective and dependent on the sensitivity of the 
measuring devices it became apparent in the generation of the spectrum in 

al 6 ' 1984? he ?i ain V°5 10 u th order recovery from SAGE, [Pisacane, et 
about 70rh At-h ld n0t be . substantiated . From intersatellite range rate, 
fl hnmi 70t K- h d t recovery is more realistic at an altitude of 400km The 
abrupt, high frequency cutoff in Figs. 11, 12, and 14 indicates there is 
information beyond 90th order in the gravity gradient components. 
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